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Abstract.  Using  a  simplified  model  of  in-orbit  radi¬ 
ance  acquisition  by  the  Advanced  Microwave  Sounding 
Unit  (AMSU-A),  we  derive  three-dimensional  temperature 
weighting  functions  for  Channel  9  measurements  (peaking 
at  ~60-90hPa)  at  all  30  cross-track  beam  positions  and 
use  them  to  investigate  the  sensitivity  of  these  radiances  to 
gravity  waves.  The  vertical  widths  of  the  weighting  func¬ 
tions  limit  detection  to  waves  with  vertical  wavelengths  of 
>10  km,  with  slightly  better  vertical  wavelength  sensitivity 
at  the  outermost  scan  angles  due  to  the  limb  effect.  Fourier 
Transforms  of  two-dimensional  cross-track  weighting  func¬ 
tions  reveal  optimal  sensitivity  to  cross-track  wavelengths  at 
the  near-nadir  scan  angles,  where  horizontal  measurement 
footprints  are  smallest.  This  sensitivity  is  greater  for  the 
AMSU-A  on  the  Aqua  satellite  than  for  the  identical  instru¬ 
ments  on  the  NOAA  meteorological  satellites,  due  to  a  lower 
orbit  altitude  and  thus  smaller  horizontal  footprints  from  an¬ 
tenna  spreading.  Small  cross-track  asymmetries  in  the  ra¬ 
diance  response  to  gravity  waves  are  found  that  peak  at  the 
mid-range  scan  angles,  with  more  symmetric  responses  at 
near-nadir  and  far  off-nadir  scan  angles.  Three-dimensional 
simulations  show  gravity  wave  oscillations  imaged  in  hori¬ 
zontal  AMSU-A  radiance  maps  swept  out  by  the  scan  pattern 
and  satellite  motion.  A  distorting  curvature  is  added  to  im¬ 
aged  wave  phase  lines  due  to  vertical  variations  in  weighting 
function  peaks  with  cross-track  scan  angle.  This  wave  dis¬ 
tortion  is  analogous  to  the  well-known  “limb  darkening”  and 
“limb  brightening”  of  microwave  radiances  acquired  from 
purely  vertical  background  temperature  profiles  by  cross¬ 
track  scanners.  Waves  propagating  along  track  are  more  vis¬ 
ible  in  these  images  at  the  outermost  scan  angles  than  those 
propagating  cross  track,  due  to  oversampling  and  narrower 
widths  of  the  horizontal  measurement  footprints  in  the  along 
track  direction.  Based  on  nominal  noise  floors  and  repre- 


Correspondence  to:  S.  D.  Eckermann 
(stephen.eckermann@nrl.navy.mil) 


sentative  lower  stratospheric  wave  temperature  amplitudes, 
our  modeling  indicates  that  Channel  9  AMSU-A  radiances 
can  resolve  and  horizontally  image  gravity  waves  with  hori¬ 
zontal  wavelengths  of  >150  km  and  vertical  wavelengths  of 
>10  km. 


1  Introduction 

While  satellite  remote  sensors  have  revolutionized  our  under¬ 
standing  of  the  global-scale  dynamics  and  chemistry  of  the 
stratosphere,  until  recently  these  instruments  have  lacked  the 
necessary  horizontal  and  vertical  resolutions  to  resolve  grav¬ 
ity  waves.  Thus,  the  observational  record  on  stratospheric 
gravity  wave  dynamics  to  date  has  relied  mostly  on  sub¬ 
orbital  observations  at  scattered  locations  around  the  globe 
(e.g.,  Eckermann  et  al.,  1995;  Wang  and  Geller,  2003).  As  a 
result,  our  understanding  of  gravity  waves  on  a  global-scale 
is  considerably  poorer  than  of  the  larger-scale  stratospheric 
dynamics  (Fritts  and  Alexander,  2003).  Since  gravity  wave 
dynamics  drive  important  aspects  of  the  global  stratospheric 
circulation,  climate  and  chemical  state  (e.g.,  Alexander  and 
Rosenlof,  2003;  Mann  et  al.,  2005),  this  lack  of  data  repre¬ 
sents  an  important  gap  in  our  knowledge. 

During  the  mid-to-late  1990s,  advances  in  limb  sound¬ 
ing  technology  at  infrared,  microwave  and  radio  wavelengths 
yielded  high-resolution  satellite  radiances  that,  for  the  first 
time,  explicitly  resolved  some  longer  wavelength  gravity 
wave  oscillations  in  the  stratosphere  (e.g.,  Fetzer  and  Gille, 
1994;  Wu  and  Waters,  1996b;  Eckermann  and  Preusse,  1999; 
Tsuda  et  al.,  2000).  While  these  data  have  provided  valu¬ 
able  glimpses  into  the  global  morphology  of  stratospheric 
gravity  waves,  they  have  also  proven  challenging  to  analyze. 
The  narrow  observational  wavelength  windows  within  which 
gravity  waves  are  visible  to  these  instruments  vary  with  chan¬ 
nel,  orbit  position  and  viewing  direction,  while  the  waves 
themselves  vary  in  wavelength  as  they  propagate  through  and 
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are  refracted  by  the  background  atmosphere.  Thus  different 
waves  are  constantly  moving  into  and  out  of  the  instrument’s 
“visibility  window,”  complicating  interpretation  of  measured 
wave  variances,  particularly  their  variations  with  time,  alti¬ 
tude  and  geographical  location  (Alexander,  1998;  Preusse  et 
ah,  2000;  Jiang  et  ah,  2004;  Wu  et  ah,  2006). 

Furthermore,  the  satellite  gravity  wave  data  acquired  to 
date  often  resemble  a  denser  global  distribution  of  subor¬ 
bital  measurements,  in  the  sense  of  providing  only  one¬ 
dimensional  cross  sections  through  three-dimensional  wave 
fields.  For  example,  some  limb  sensors  return  high- 
resolution  vertical  temperature  profiles  with  wave  oscilla¬ 
tions  superimposed,  but  separated  too  far  apart  horizon¬ 
tally  to  resolve  horizontal  structure  (Eckermann  and  Preusse, 
1999;  Tsuda  et  ah,  2000).  Thus  these  data  provide  only  ver¬ 
tical  profiles  of  wave  oscillations,  similar  to  radiosonde  and 
rocketsonde  data.  Similarly,  limb-tracking  measurements 
from  the  Microwave  Limb  Sounder  (MLS)  on  the  Upper  At¬ 
mosphere  Research  Satellite  (UARS)  return  high-resolution 
radiances  from  a  given  altitude  that  resolve  wave  fluctua¬ 
tions  along  the  orbital  track  (Wu  and  Waters,  1996b;  Jiang 
et  ah,  2004),  similar  to  in  situ  aircraft  data.  Many  funda¬ 
mental  gravity  wave  properties,  such  as  vertical  fluxes  of 
horizontal  pseudomomentum  densities,  cannot  be  quantified 
from  one-dimensional  profiles  alone.  Thus,  we  require  new 
satellite  observations  that  can  resolve  gravity  wave  oscilla¬ 
tions  in  two  or  three  spatial  dimensions.  Here  we  investigate 
whether  two-dimensional  maps  of  lower  stratospheric  radi¬ 
ances  acquired  by  the  Advanced  Microwave  Sounding  Unit- 
A  (AMSU-A)  can  horizontally  image  gravity  waves. 

AMSU-A,  described  in  Sect.  2,  is  the  latest  in  a  series  of 
cross-track  scanning  satellite  microwave  radiometers.  It  con¬ 
solidates  over  30  years  of  accumulated  experience  and  tech¬ 
nological  development  in  passive  microwave  remote  sens¬ 
ing  of  atmospheric  temperatures  from  space,  beginning  with 
the  Nimbus  E  Microwave  Sounder  (NEMS)  on  the  Nim¬ 
bus  5  satellite.  NEMS  included  three  narrowband  tempera¬ 
ture  channels  within  the  50-60  GHz  O2  thermal  band  whose 
vertical  weighting  functions  peaked  at  ~4,  10  and  17  km  al¬ 
titude  (Waters  et  ah,  1975).  While  NEMS  viewed  only  in 
the  nadir,  it  was  realized  that  the  receiving  antenna  could  be 
scanned  across  the  satellite  track  to  acquire  radiances  at  a 
series  of  other  atmospheric  locations  either  side  of  the  sub¬ 
satellite  point.  The  first  instrument  with  this  scanning  capa¬ 
bility  was  the  Scanning  Microwave  Spectrometer  (SCAMS) 
on  Nimbus  6,  which  performed  sequential  “step  and  stare” 
measurements  at  13  cross-track  off-nadir  viewing  angles  dis¬ 
tributed  symmetrically  out  to  ±1200  km  either  side  of  the 
subsatellite  point  (Grody  and  Pellegrino,  1977).  SCAMS  was 
superseded  by  the  Microwave  Sounding  Unit  (MSU),  which 
acquired  data  in  a  similar  fashion  at  11  cross-track  beam 
positions.  MSU  instruments  flew  on  the  NOAA-6  through 
NOAA-14  meteorological  satellites. 

SCAMS  and  MSU  used  antennas  with  beamwidths  that 
translated  into  half-power  horizontal  measurement  footprints 


at  the  ground  with  diameters  of  >100  km  at  nadir  and  nearer 
200-300  km  at  the  off-nadir  angles  (Grody  and  Pellegrino, 
1977;  Kidder  and  Vonder  Haar,  1995).  Since  extratropical 
gravity  waves  have  horizontal  wavelengths  in  the  ~5-500  km 
range,  these  instruments  could  not  horizontally  image  grav¬ 
ity  wave  perturbations  in  their  acquired  radiances,  since  their 
footprint  “pixels”  were  too  broad.  Conversely,  the  AMSU-A 
receiving  antenna  has  a  narrower  antenna  polar  diagram  that 
yields  smaller  surface  footprints  (see  Fig.  1  of  Kidder  et  ah, 
2000).  Furthermore,  it  has  many  more  channels,  including  6 
stratospheric  temperature  channels,  and  improved  radiomet¬ 
ric  accuracy.  All  these  improvements  over  MSU  mean  that 
AMSU-A  might  resolve  and  image  long  wavelength  gravity 
waves  oscillations  in  its  stratospheric  radiances. 

Wu  (2004)  investigated  this  possibility  experimentally 
by  extracting  and  analyzing  along-track  perturbations  in 
AMSU-A  stratospheric  radiances  at  various  cross-track  scan 
angles.  Maps  of  along-track  Channel  13  radiance  variances 
(peaking  at  ~5  hPa)  in  the  Southern  Hemisphere  showed  en¬ 
hancements  over  land  masses  and  within  strong  vortex  winds 
that  agreed  well  with  variances  in  along-track  UARS  MLS 
limb  radiances  acquired  at  a  similar  altitude  and  time  of  year. 
This  strongly  suggests  that  AMSU-A  stratospheric  radiances 
do  resolve  some  gravity  wave  oscillations. 

Wu  and  Zhang  (2004)  used  this  analysis  technique  to 
isolate  AMSU-A  stratospheric  radiance  perturbations  dur¬ 
ing  19-21  January  2003,  when  strong  baroclinic  deforma¬ 
tion  of  tropospheric  jet  streams  near  the  east  coast  of  the 
United  States  appeared  to  radiate  inertia  gravity  waves  into 
the  stratosphere.  On  plotting  radiance  perturbations  in  all  30 
beams  as  horizontal  “pushbroom”  images,  two-dimensional 
wavelike  oscillations  were  revealed  with  impressively  coher¬ 
ent  linear  phase  lines  in  the  horizontal.  These  structures  were 
seen  in  different  AMSU-A  stratospheric  channels,  extend¬ 
ing  from  the  lower  to  the  upper  stratosphere.  Their  results 
appear  to  show  AMSU-A  “imaging”  the  horizontal  struc¬ 
ture  of  a  jet-generated  gravity  wave  packet  as  it  propagates 
through  the  stratosphere.  Such  measurements,  if  validated, 
could  provide  a  valuable  new  global  measurement  capabil¬ 
ity  for  stratospheric  gravity  waves,  since  to  date  only  some 
very  limited  4.3  /xm  Midcourse  Space  Experiment  (MSX)  ra¬ 
diances  have  horizontally  imaged  stratospheric  gravity  waves 
from  orbit  (Dewan  et  ah,  1998). 

Experience  to  date  has  shown  that  a  detailed  understand¬ 
ing  of  the  visibility  characteristics  of  a  satellite  instrument  to 
gravity  waves  is  critical  to  a  proper  interpretation  of  wave- 
induced  perturbations  in  the  acquired  data  (Alexander,  1998; 
McLandress  et  ah,  2000;  Preusse  et  ah,  2000;  Jiang  et  ah, 
2004;  Wu  et  ah,  2006).  How  gravity  waves  might  manifest 
in  stratospheric  radiances  acquired  by  AMSU-A  is  not  well 
understood  at  present.  Thus,  in  Sect.  3  we  develop  a  simple 
forward  model  of  the  in-orbit  AMSU-A  Channel  9  radiance 
acquisition,  aimed  at  understanding  how  three-dimensional 
gravity  wave  oscillations  might  appear  in  these  radiances. 
While  we  take  care  to  validate  the  accuracy  of  the  model 
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we  develop,  we  also  seek  relevant  simplifications  whereever 
possible  so  that  we  can  gain  physical  insight  into  how  grav¬ 
ity  waves  are  observed,  and  can  also  apply  the  model  more 
easily  to  interpret  future  AMSU-A  observations.  We  quan¬ 
tify  the  sensitivity  of  Channel  9  radiances  to  gravity  waves 
of  different  vertical  and  horizontal  wavelengths  via  spectral 
analysis  of  the  resulting  weighting  functions  (Sect.  4)  and 
three-dimensional  forward  modeling  using  idealized  three- 
dimensional  gravity  wave  temperature  oscillations  (Sect.  5). 
An  observational  study  of  a  gravity  wave  event  imaged  in 
AMSU-A  radiances  that  utilizes  and  validates  the  modeling 
predictions  from  this  work  forms  the  subject  of  a  companion 
paper  (Eckermann  et  al.,  2006b). 


2  Advanced  Microwave  Sounding  Unit-A 

The  Advanced  Microwave  Sounding  Unit  (AMSU)  is  a 
cross-scanning  passive  microwave  sounding  instmment  cur¬ 
rently  deployed  on  five  different  satellites:  the  NOAA-15, 
NOAA-16,  NOAA-17  and  NOAA-18  meteorological  satel¬ 
lites  (Mo,  1996;  Kidder  et  ah,  2000),  and  NASA’s  Earth  Ob¬ 
servation  System  (EOS)  Aqua  satellite  (Lambrigtsen,  2003). 

AMSU  consists  of  three  separate  hardware  modules,  de¬ 
noted  AMSU-A1,  AMSU-A2,  and  AMSU-B.  AMSU-B  is 
known  as  the  Humidity  Sounder  for  Brazil  (HSB)  on  EOS 
Aqua,  and  was  superseded  by  the  Microwave  Humidity 
Sounder  (MHS)  on  NOAA-18:  we  do  not  use  data  from  these 
modules  in  this  study.  AMSU-A1  has  12  channels  sampling 
thermal  oxygen  emissions  in  the  50-58  GHz  band,  and  one 
channel  at  89  GHz,  while  AMSU-A2  has  two  water  vapor 
channels  at  23.8  and  3 1.4  GHz:  together,  these  two  mod¬ 
ules  comprise  AMSU-A.  The  two  AMSU-A2  channels  are 
denoted  AMSU-A  channels  1  and  2.  AMSU-A  channels 
3-8  are  the  first  six  AMSU-A1  channels  centered  at  50.3, 
52.8,  53.596±0. 1 15,  54.4,  54.94  and  55.5  GHz,  respectively, 
with  channel  5  sampling  two  O2  wing  emission  lines  either 
side  of  the  central  line.  Channels  9-14  sample  O2  wing  line 
emissions  centered  at  57.290  GHz,  and  channel  15  samples 
emissions  at  89  GHz.  AMSU-A  channels  1-3  and  15  are 
“window”  channels  that  mostly  sense  surface  parameters. 
The  other  11  channels  (4—14)  are  atmospheric  temperature 
sounding  channels  that  acquire  radiances  from  progressively 
higher  in  the  atmosphere,  from  ~900  hPa  for  channel  4,  up 
to  ~2.5  hPa  for  channel  14.  Channels  4-8  peak  in  the  tro¬ 
posphere,  while  channels  9-14  are  stratospheric.  For  further 
details,  see  Table  1  and  Fig.  3  of  Goldberg  et  al.  (2001). 

Figure  1  depicts  the  AMSU-A  scanning  pattern  with  re¬ 
spect  to  the  orbital  geometry.  The  satellite,  depicted  by 
the  square  on  the  Z  axis  at  its  orbit  altitude  Zsat,  or¬ 
bits  along  the  X  axis,  tracing  out  a  ground  track  at  T=0. 
The  NOAA  satellites  orbit  the  Earth  every  102  min  at  an 
altitude  Zsat=833km,  and  so  move  along  the  X  axis  at 
+7.4kms_1.  EOS  Aqua  orbits  at  a  somewhat  lower  alti¬ 


Fig.  1.  Schematic  depiction  of  the  AMSU-A  measurements.  The 
satellite  is  depicted  by  the  square  on  the  Z  axis  and  orbits  along  the 
X  axis.  From  this  platform.  AMSU-A  performs  j  =  1. .  .30  sequen¬ 
tial  step-and-stare  measurements  at  equispaced  off-nadir  cross-track 
scan  angles  fi  j,  ranging  from  ji\  =— 48.33°  to  ^3o=+48.33°.  Sign 
convention  for  quantities  on  the  Y  axis  was  chosen  to  preserve  a 
right-handed  coordinate  system  (X,  Y,  Z). 


tude  of  Zsat=705km  and  has  an  orbital  period  of  98.8  min 
(Parkinson,  2003),  yielding  a  7.5  km  s-1  satellite  velocity. 

AMSU-A  performs  periodic  cross-track  scans  along  the 
Y  axis.  As  depicted  in  Fig.  1,  each  scan  cycle  consists  of 
30  individual  “step  and  stare”  measurements  at  30  different 
beam  positions  j  defined  by  their  different  equispaced  off- 
nadir  cross-track  viewing  angles  ft  j  with  respect  to  the  sub¬ 
satellite  point,  given  by 


-155  +  10; 


deg.  0  =  1...  30). 


(1) 


Thus  each  scan  starts  at  beam  position  ;  =  1 
(/li=— 48.33°),  then  progresses  along  Y  in  3.33°  steps  to  a 
final  measurement  at  beam  position  j  =  30  (/l3o=+48.33°): 
see  also  Fig.  2.  The  scanning  pattern  is  symmetric  about  the 
subsatellite  ground  track  (Y  —0).  Each  staring  measurement 
lasts  for  0.165  s  followed  by  ~37.5  ms  to  rotate  to  the  next 
beam  position,  yielding  a  measurement  every  Ml.  2025  s 
and  ~6  s  to  perform  a  complete  cross-track  scan  of  30 
measurements.  After  each  such  scan  cycle,  ~2  s  is  devoted 
to  radiometric  recalibration  via  viewing  of  an  internal 
blackbody  source  and  cold  space,  so  that  the  total  scan  cycle 
period  is  8  s.  For  further  details,  see  Lambrigtsen  (2003). 


3  Simple  forward  model  of  AMSU-A  Channel  9  radi¬ 
ance  acquisition 

3.1  Absorption 

Wu  (2004)  presented  a  simple  model  of  the  AMSU-A  radi¬ 
ance  measurement,  which  we  develop  further  here.  We  fo¬ 
cus  here  on  the  lowest-altitude  stratospheric  channel  (Chan¬ 
nel  9),  which  measures  the  central  57.290  GHz  O2  wing  line 
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beam  position  j 


Fig.  2.  Asterisks  connected  with  solid  line  show  the  off-nadir  cross¬ 
track  scan  angles  fl  j  of  AMSU-A  beam  positions  7  =  1-15  depicted 
in  Fig.  1.  Diamonds  connected  with  dotted  line  show  the  corre¬ 
sponding  off-nadir  angles  of  the  line-of-sight  ray  to  the  satellite 
at  the  “measurement  point,”  defined  as  the  point  of  peak  response 
in  the  3-D  weighting  functions  (see  later).  These  values  are  derived 
from  Eq.  (A2)  in  Appendix  A  for  a  nominal  measurement  altitude  of 
Zq= 18km.  These  p:  angles  differ  slightly  from  fi;  at  the  farthest 
off-nadir  scan  angles  due  to  Earth  curvature:  see  Eq.  (A2). 


emission.  Initially,  we  will  assume  that  this  emission  peaks 
low  enough  in  the  atmosphere  (~60-90  hPa)  that  pressure¬ 
broadening  dominates,  so  that  it  can  be  modeled  to  a  good 
approximation  as  a  single  Lorentz  line.  Then,  the  absorption 
coefficient  is  the  Lorentz-line  solution  for  a  gas  with  constant 
mixing  ratio:  kv=Bp,  where  p  is  atmospheric  pressure  and 
B  is  a  line  constant  (Houghton  et  ah,  1984).  Radiance  ab¬ 
sorption  along  a  ray  path  distance  s  from  the  Earth’s  surface 
is  then  given  by 


r(s)  = 


[  Bp{s')po2{s')ds'  =  f 
Jo  Jo 

Ap2{s') 


-Bq02P2(s') 


o  gH(s') 


ds' 


_  r  Air(S')  ,  r* 
Jo  H(s')  Jo 


dr  {s') 
ds 


ds  , 


(2) 


where  po2  is  the  O2  mass  density,  <702  is  the  (constant)  O2 
mass  mixing  ratio,  and  H  is  pressure  scale  height  (Houghton 
et  ah,  1984). 


3.2  Antenna  pattern 


Since  the  central  AMSU-A  antenna  beam  has  an  approxi¬ 
mately  Gaussian-shaped  polar  diagram  (Lambrigtsen,  2003), 
we  approximate  the  polar  diagram  along  the  cross  track  scan 
axis  Y  at  beam  position  j  as 


Pi-d(P,  Pj)  oc  exp 


Pw 


Phpbw 
2(ln2)1/2  ’ 


(3) 

(4) 
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where  p  is  the  off-nadir  angle  along  the  cross-track  Y  direc¬ 
tion,  and  the  Gaussian  antenna  width  Pw  is  related  via  Eq.  (4) 
to  the  quoted  full-width  half-power  (3  dB)  beamwidth  of  the 
AMSU-A1  antenna  /3hpbw-  which  is  nominally  3.3°  but  was 
measured  by  Mo  (1996)  to  be  ~3.51°  for  Channel  9.  We  use 
this  latter  value  in  our  model.  Since  the  antenna  pattern  is 
symmetric  in  the  along-track  and  cross-track  directions  (Mo, 
1999),  our  two-dimensional  antenna  pattern  is 

P2-d(«,  P,  Pj)  cx  Pi-D(P,  Pj)Pi-D(a ,  0),  (5) 

where  a  is  the  off-nadir  angle  in  the  along-track  X  direc¬ 
tion:  note  that  aj= 0  for  all  j  in  Eq.  (5)  since  AMSU-A  scans 
cross-track  only  (see  Fig.  1). 

3.3  Radiative  transfer  model 


We  adopt  a  simplified  model  of  microwave  radiative  transfer 
in  the  ( X .  Y.  Z )  plane  that  ignores  noise,  cloud  scattering  and 
surface  effects,  and  invokes  the  Rayleigh-Jeans  approxima¬ 
tion  to  the  microwave  Planck  function  (e.g.,  Staelin,  1977). 
The  brightness  temperature  (radiance)  contribution  for  beam 
position  j  at  a  distance  s  from  the  surface  along  a  ray  path 
directed  to  the  satellite  is  quantified  by  the  gradient  term 


dTn  r+*/2  r+n/2 

—P(s,Pj)(x  /  /  P2.D(a,p,Pj ) 

ds  J-n/2  J-j t/2 

T  (. Y Y',  Z')  e~r(s)dadp , 


dr  {s) 
ds 


(6) 


where  T {X,  T,  Z)  is  the  atmospheric  temperature  and 
X'(s,  a,  P),  Y'(s,  a,  /6),  Z'(i,  a,  P)  are  the  positions  in  the 
atmosphere  in  the  ( X .  Y.  Z)  frame  in  Fig.  1  implied  by  geom¬ 
etry  for  a  line-of-sight  ray  to  the  satellite  inclined  at  off-nadir 
angles  a  and  p.  The  total  brightness  temperature  acquired  by 
the  instrument  at  beam  position  j  is 

rss*  dTn 

Tb(Pj)=  -PL(s\pj)ds'.  (7) 

Jo  ds 

where  ssat  is  the  total  length  of  this  ray  path  from  the  ground 
to  the  satellite. 

Our  model  solves  this  simplified  radiative  transfer  numer¬ 
ically  in  the  Y-Z  plane  by  integrating  along  successive  ray 
paths  s  for  a  full  range  of  angles  p  centered  about  the  cen¬ 
tral  antenna  beam  Pj  for  a—0.  These  calculations  use  a 
full  spherical  geometry  outlined  briefly  in  Appendix  A:  see 
Fig.  16. 


3.4  Weighting  functions 


For  a  purely  vertical  mean  temperature  profile  T(Z), 
our  numerical  model  calculations  in  Sect.  3.3  yield  the 
two-dimensional  (2-D)  temperature  weighting  functions 
Wj(Y,Z)  and  then  the  one-dimensional  (1-D)  vertical 
weighting  functions  W j(Z)=  ){Y,  Z)dY  for  each 

beam  position  j .  Since  our  radiative  transfer  model  is  ideal¬ 
ized,  to  simplify  the  calculations  further  we  used  H=  7.5  km 
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in  Eq.  (2),  then  adjusted  the  absorption  constant  A  in  Eq.  (2) 
to  give  best  agreement  with  the  Channel  9  1-D  vertical 
weighting  functions  for  the  near-nadir  beam  positions  7=15- 
16  published  by  Goldberg  et  al.  (2001),  based  on  rigorous 
radiative  transfer  calculations.  We  used  the  1976  U.S.  Stan¬ 
dard  Atmosphere  temperature  profile  for  T (Z),  although  un¬ 
der  these  model  approximations  our  final  weighting  func¬ 
tions  have  no  explicit  dependence  on  the  temperature  profile 
used. 

3.4.1  1-D  vertical  weighting  functions 

The  lower  dotted  curve  in  Fig.  3  shows  the  vertical  weight¬ 
ing  function  Wj  (Z)  for  the  near-nadir  beams  (7=15-16)  from 
our  simple  Lorentz-line  model  with  a  fitted  constant  A.  It  re¬ 
produces  the  Goldberg  et  al.  (2001)  profile  (gray  solid  curve) 
quite  well  at  altitudes  near  the  peak  response.  In  fact  these 
near-nadir  model  weighting  functions  have  shapes  similar  to 
the  idealized  analytical  profile  solutions  ocp2  exp[— Ap2]  ap¬ 
propriate  for  a  single  Lorentz  line  emission  from  a  gas  with 
constant  mixing  ratio  viewed  in  the  nadir  without  antenna 
spreading  (Houghton  et  al.,  1984;  Grody,  1993),  where  A  is 
a  constant  related  to  A. 

This  same  constant  A  value  is  then  used  in  our  model  to 
evaluate  all  the  other  1-D  and  2-D  weighting  functions  at  all 
beam  positions  7=1  ...  30.  The  second  (upper)  dotted  curve 
in  Fig.  3  shows  the  resulting  model-generated  1-D  vertical 
weighting  function  at  the  largest  off-nadir  beam  positions 
7=1  and  7=30.  The  agreement  with  the  corresponding  pro¬ 
file  of  Goldberg  et  al.  (2001)  is  again  quite  good  given  the 
simplicitly  of  our  model. 

Nonetheless,  the  model  results  diverge  from  the  Gold¬ 
berg  et  al.  (2001)  profiles  most  noticeably  at  the  tails  of 
the  weighting  functions  at  altitudes  away  from  the  peak  re¬ 
sponse.  Though  these  differences  seem  minor  given  the 
small  values  of  the  weighting  functions  here,  simple  1- 
D  modeling  calculations  using  sinusoidal  vertical  temper¬ 
ature  perturbations  revealed  significant  differences  in  ac¬ 
quired  brightness  temperature  perturbation  amplitudes  be¬ 
tween  these  profiles  and  the  Goldberg  et  al.  (2001)  profiles 
for  shorter  vertical  wavelength  perturbations. 

To  obtain  better  fits,  we  tuned  our  absorption  coefficient 
A  in  Eq.  (2)  to  have  a  pressure-altitude  dependence  A(Z) 
that  parameterizes  the  net  effects  of  non-Lorentzian  and/or 
overlapping  lines  (following,  e.g..  Waters  et  al.,  1975;  Poon, 
1980;  Grody,  1993).  Closer  fits  were  obtained  by  reducing 
A  with  pressure  altitude  approximately  linearly  below  20  km 
and  above  25  km.  Solid  curves  in  Fig.  3  show  the  verti¬ 
cal  weighting  functions  produced  by  these  “tuned”  Lorentz 
model  runs.  Since  they  fit  the  Goldberg  et  al.  (2001)  profiles 
better  at  all  altitudes,  we  use  weighting  functions  from  this 
tuned  model  in  all  subsequent  analysis. 

The  off-nadir  W j(Z)  curves  in  Fig.  3  peak  higher  in  alti¬ 
tude  since  they  traverse  longer  pathlengths  through  the  atmo¬ 
sphere:  the  so-called  “limb  effect”  (e.g.,  Grody,  1993).  Note 


NO  A  A  15-18  AMSU-A  Channel  9 


normalized  response 


Fig.  3.  Channel  9  1-D  vertical  temperature  weighting  functions 
Wj(Z)  for  the  smallest  off-nadir  angle  of  ^=±1.67°  and  largest 
off-nadir  angle  of  ^=±48.33°,  derived  from  our  simple  Lorentz 
model  (constant  A;  dotted  cuve)  and  tuned  Lorentz  model  (height- 
varying  A,  solid  curve).  Gray  curves  show  results  of  Goldberg  et  al. 
(2001)  using  a  complete  radiative  transfer  model. 

also  that  the  off-nadir  Channel  9  weighting  function  is  nar¬ 
rower  vertically  and  thus  has  a  larger  maximum  value.  This 
indicates  that,  for  infinitely  long  horizontal  wavelengths,  off- 
nadir  beams  should  be  slightly  more  sensitive  to  vertical  tem¬ 
perature  oscillations  than  the  near-nadir  beams. 

3.4.2  2-D  and  3-D  weighting  functions 

Figure  4  profiles  Wj(Y,  Z)  from  our  tuned  Lorentz  model 
for  AMSU-A  measurements  from  the  NOAA  satellites  at  all 
15  different  off-nadir  scan  angles.  The  results  illustrate  not 
only  an  increase  in  the  peak  altitude  of  the  weighting  func¬ 
tion  with  increasing  \fij  due  to  the  limb  effect,  but  also  an 
increasing  cross-track  width.  These  features  are  studied  fur¬ 
ther  in  Sect.  3.5. 

Since  AMSU-A  does  not  scan  along-track,  we  assume  no 
significant  limb  effect  along-track,  so  that  the  along-track 
variations  in  the  weighting  functions  are  controlled  solely 
by  the  antenna  pattern.  This  yields  a  separable  form  for  our 
assumed  three-dimensional  (3-D)  AMSU-A  weighting  func¬ 
tion 

Wj  (X,Y,Z)  =  G(X)Wj(Y,Z),  (8) 

where  G(X)  specifies  the  along-track  spreading  at  a  given 
T,  Z  and  beam  position  j  controlled  by  the  antenna  pattern 
Eq.  (5)  and  spherical  geometry.  There  will  be  a  small  ad¬ 
ditional  smearing  contribution  along  the  X  direction  due  to 
the  7.4 kms-1  motion  of  the  satellite  and  the  0.165  s  staring 
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Fig.  4.  Modeled  AMSU-A  Channel  9  weighting  functions 
Wj(Y,  Z)  for  beam  positions  j= 1 ...  15.  Thick  contours  show  the 
half-power  (50%)  levels,  the  outer  ranges  of  which  are  projected 
to  the  surface  to  depict  the  horizontal  footprints.  Other  contours 
show  the  70%,  90%  and  99%  levels.  Dotted  lines  show  line  of  sight 
ray  paths  from  the  peaks  in  each  weighting  function  to  the  satellite 
point,  taking  into  account  Earth  curvature:  the  off  nadir  angles 
that  each  of  these  ray  paths  makes  at  the  weighting  function  peak 
are  plotted  as  the  dotted  curve  in  Fig.  2.  The  W:(Y,  Z)  contours  are 
plotted  along  a  vertical  plane  aligned  along  an  approximate  scan 
axis  that  is  offset  slightly  from  the  cross-track  axis  Y  due  to  the 
7.4  km  s_1  motion  of  the  satellite  along  the  X  axis  during  the  3  s 
taken  to  scan  through  beam  positions  j=l~15.  The  actual  X-Y  scan 
track  is  plotted  in  green,  and  curves  slightly  in  the  X-Y  plane.  Note 
the  different  X  and  Y  axis  scales:  footprint  diameters  along  track 
are  in  fact  narrower  than  those  cross  track  (see  Fig.  5). 


time,  but  the  effect  is  small  in  comparison  to  G(X)  and  so  is 
ignored  hereafter. 

3.5  Horizontal  footprints 


The  half-power  levels  for  Wj (Y,  Z)  (thick  contours  in  Fig.  4) 
are  used  to  specify  the  cross-track  diameters  of  the  mea¬ 
surement  “footprints”  (fy)j-  Using  Eq.  (8),  and  given  the 
symmetric  antenna  pattern  of  Eqs.  (3)-(4),  we  can  convert 
these  (fy)  j  values  from  our  model  into  associated  along- 
track  half-power  footprint  widths  (fx)j  via  approximate  ge¬ 
ometrical  relations  linking  the  two  for  a  curved  Earth:  see 
Appendix  A.  The  result  is  a  scaling  ratio 


(fx)j 
( fr)j 


2  sin  4>(fi  j)  tan 

/Shpbw 

L  2  J 

sin  Pj 

[m  +  te2sw) 

<H  Pj  - 

Aipbw  i 

2  ’ 

(9) 


Fig.  5.  Ratio  of  along-track  to  cross-track  footprint  diameters 
( fx) j Hfy ) j  versus  scan  angle  fij  based  on  Eqs.  (9)  and  (10)  for 
NOAA  satellite  orbit  parameters  and  a  nominal  Channel  9  measure¬ 
ment  altitude  of  Zc=18  km. 


where 


(p(/3)  =  arcsin 


(Re  +  Zgat)  sin/1 

Re  +  Zc 


(10) 


is  the  angle  from  the  center  of  the  Earth  between  the  satellite 
and  the  measurement  point  for  a  given  cross-track  off-nadir 
angle  f}.  Re  is  Earth  radius  and  Zq  is  the  altitude  of  the 
measurement  point.  Given  (fy)  j  from  our  model,  we  de¬ 
rive  (fx)j  using  Eqs.  (9)  and  (10).  Ratios  from  Eq.  (9)  are 
plotted  in  Fig.  5  for  NOAA  satellite  orbit  parameters.  The 
elliptical  half-power  horizontal  footprints  implied  by  (fx)j 
and  (  fy)  j  are  projected  onto  the  Z=0  surface  in  Fig.  4.  Fig¬ 
ure  6a  shows  how  the  7.4  km  s~ 1  velocity  of  the  NOAA  satel¬ 
lites  along  the  X  axis  maps  these  observational  footprints  into 
two-dimensional  horizontal  measurement  coverage  as  the  in¬ 
strument  cyclically  scans  cross-track.  Almost  all  regions  at 
cross-track  distances  ±  1100  km  either  side  of  the  satellite 
ground  track  are  measured  with  >50%  sensitivity.  Footprint 
sizes  increase  at  the  outer  scan  angles  and  overlap  (oversam¬ 
ple)  in  the  X  direction.  These  footprints,  derived  from  our 
tuned  Channel  9  radiance  acquisition  model,  compare  well 
with  those  from  more  detailed  calculations  of  surface  foot¬ 
prints  for  AMSU-A  window  channels  (see,  e.g..  Fig.  1  of 
Kidder  et  al.,  2000). 

Figure  6b  shows  the  corresponding  footprints  calculated 
for  the  AMSU-A  on  EOS  Aqua.  The  lower  orbit  altitude 
(Zsat= 705  km  compared  to  Zsat=833  km  for  the  NOAA  satel¬ 
lites)  produces  correspondingly  smaller  footprints  and  less 
total  horizontal  coverage  cross  track.  Since  the  satellite  ve¬ 
locity  is  about  the  same,  these  smaller  footprints  mean  that 
not  all  regions  within  the  swath  are  sampled  at  >50%  sensi¬ 
tivity  from  EOS  Aqua. 
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(a)  NOAA  15-18  Channel  9  AMSU-A  Footprints 


cross-track  distance  Y  (km) 


(b)  EOS  Aqua  Channel  9  AMSU-A  Footprints 


cross-track  distance  Y  (km) 


Fig.  6.  AMSU-A  Channel  9  horizontal  footprints  as  a  function  of 
along-track  and  cross-track  distances  traced  out  by  the  AMSU-A 
scanning  pattern  from  (a)  NOAA-15  through  NOAA-18  satellites 
and  (b)  EOS  Aqua.  Orange  line  shows  the  satellite  ground  track. 
Green  curves  show  the  scanning  pattern  from  right-to-left  across  Y 
as  the  satellite  moves  along  X. 


3.6  Simulated  brightness  temperature  measurement 

Having  derived  weighting  functions,  we  now  use  them  to 
simulate  radiance  acquisition  in  the  presence  of  varying  at¬ 
mospheric  temperatures.  If  we  center  the  peak  sensitivity  in 
the  various  3-D  weighting  functions  Wj(X,  Y,  Z)  at  X—Y— 0 
and  define  X j,  Y /  to  be  the  current  central  position  of  this 
beam’s  horizontal  measurement  footprint,  then  the  acquired 
brightness  temperature,  Eq.  (7),  can  be  reexpressed  as 

/oo  poo  pZsat 

/  /  Wj(X-Xj,Y-Yj,Z ) 

-oo  J  —oo  J 0 

T(X,  Y,  Z)dXdYdZ,  (11) 

where  T (X,  Y,  Z)  is  atmospheric  temperature  and  here  the 
3-D  weighting  functions  are  normalized  such  that 

/OO  /’OO  p  Zsat 

/  /  Wj(X,  Y,  Z)dXdYdZ  —  1.  (12) 

-oo  J  —oo  J0 


We  evaluate  Eq.  (11)  numerically  in  Fig.  7  in  simulat¬ 
ing  an  AMSU-A  measurement  from  NOAA  satellites  of  the 
mean  vertical  temperature  profile  over  southern  Scandinavia 
on  14  January  2003  at  12:00  UTC,  as  specified  by  analysis 
fields  in  Fig.  7a  from  the  European  Center  for  Medium  Range 
Weather  Forecasts  (ECMWF).  Figure  7b  plots  the  simulated 
Channel  9  brightness  temperatures,  and  Fig.  7c  plots  the 
mean  brightness  temperatures  as  a  function  of  scan  angle  fij. 

We  see  a  decrease  in  the  Channel  9  brightness  temperature 
with  increasing  caused  by  the  increase  in  height  of  the 
peak  sensitivity  of  the  weighting  functions  with  increasing 
l/l/l  (Figs.  3^1)  coupled  with  the  decrease  with  height  in  tem¬ 
peratures  in  Fig.  7a  in  the  ~60-90hPa  region  where  Chan¬ 
nel  9  measurement  sensitivities  peak  (Fig.  3).  This  simulates 
the  well-known  “limb  darkening”  effect  typically  observed 
in  raw  swath-scanned  microwave  radiance  imagery  from  the 
troposphere  due  to  negative  tropospheric  lapse  rates  (Kidder 
and  Vonder  Haar,  1995):  such  effects  are  corrected/adjusted 
when  retrieving  temperatures  from  these  data  (Goldberg  et 
al.,  2001).  Figure  7c  reveals  a  fairly  smooth  cross-track  trend 
in  the  limb  darkening  with  increasing  1/3/1 ,  consistent  with 
the  observational  study  of  Wu  (2004),  who  found  that  cross¬ 
track  trends  in  raw  AMSU-A  brightness  temperatures  due  to 
limb  effects  could  be  accurately  removed  by  fitting  a  least- 
squares  polynomial  to  the  cross-track  radiances. 


4  Spectral  sensitivity  to  gravity  waves 

From  Eq.  (1 1)  we  see  that  (to  within  a  sign  and  vertical  shift 
convention)  the  3-D  AMSU-A  weighting  functions  are  con¬ 
volved  through  3-D  atmospheric  temperature  fields  by  the 
cyclical  scanning  pattern  to  yield  a  final  horizontal  map  of 
brightness  temperatures.  Since  convolution  corresponds  to 
multiplication  in  the  Fourier  domain,  then  the  Fourier  Trans¬ 
forms  of  the  weighting  functions  Wj(X,Y,  Z)  define  the 
spectral  visibility  (Jiang  et  al.,  2004)  of  AMSU-A  radiances 
to  sinusoidal  temperature  structure  with  different  horizontal 
and  vertical  wavelengths. 

4.1  Fourier  Transforms  of  two-dimensional  weighting 
functions 

Figure  8  plots  the  Fourier  Transforms  W j(ky ,  kz)  of  the 
NOAA  15-18  AMSU-A  weighting  functions  Wj(Y,  Z)  at 
a  sequence  of  beam  positions,  ranging  from  far  off-nadir 
(panel  a,  j= 2,  Pj=— 45°)  to  near-nadir  (panel  d,  j=  1 5, 
Pj  =  — 1.67°).  For  a  given  gravity  wave  of  wavenumber 
{ky ,  kz).  these  2-D  spectra  define  visibilities  at  both  the 
beam  position  j  and  its  conjugate  beam  position  j*=3 1  —  j 
(such  that  Pjt=—/3j),  since 

Wjt(ky,  kz)  =  Wj (  ky ,  kz).  (13) 

Equation  (13)  arises  due  to  the  symmetry  of  the  AMSU-A 
scanning  pattern  about  the  subsatellite  point  (ignoring  small 
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Fig.  7.  (a)  Mean  temperature  profile  from  the  ECMWF  analysis  for  14  January  2003  at  12:00  UTC  at  10°  E,  60°  N;  (b)  model  simulation 
of  the  AMSU-A  Channel  9  brightness  temperatures  Tg{Xj,  Y  j)  obtained  from  sampling  this  temperature  profile  using  Eq.  (11);  (c)  mean 
brightness  temperatures  from  (b)  as  a  function  of  f)j. 


asymmetries  introduced  by  the  satellite  motion:  see  Fig.  6). 
Physically,  Eq.  (13)  means  that  observing  an  infinite  wave 
train  of  wavenumber  (ky ,  kz)  at  beam  position  7*  is  equiva¬ 
lent  to  observing  the  identical  antipropagating  wave  train  of 
wavenumber  (— ky ,  kz)  at  beam  position  j .  Note  also  that 
W j(ky,  kz)—Wj(—ky ,  —  kz)  for  all  ky,  kz  and  j . 

Since  the  near-nadir  observations  have  much  narrower 
cross-track  footprint  diameters  (fy)j  than  the  far  off-nadir 
observations,  they  can  resolve  structure  with  much  shorter 
cross-track  wavelengths  Xy=2n /\ky  |,  as  evidenced  by  the 
broader  spectral  sensitivity  along  the  ky  axis  in  Fig.  8d  com¬ 
pared  to  Fig.  8a.  The  sensitivity  to  vertical  wavelengths 
Xz—2n /\kz\  is  slightly  better  for  the  /'= 2  spectrum  com¬ 
pared  to  the  7=15  spectrum  due  to  the  vertical  narrowing  of 
the  weighting  function  response  due  to  the  limb  effect,  noted 
earlier  in  Fig.  3. 

Gravity  wave  detection  requires  a  brightness  temperature 
oscillation  that  lies  above  the  noise  floor,  which  is  quantified 
for  AMSU-A  by  the  so-called  noise  equivalent  delta  temper¬ 
ature  (NEAT).  An  NEAT=0.236K  was  measured  prior  to 
launch  for  Channel  9  of  the  NOAA  AMSU-A  instruments 
(Mo,  1996).  However,  Wu  (2004)  used  observed  minima  in 
AMSU-A  brightness  temperature  variances  on  NOAA  15,  16 
and  17  to  infer  a  smaller  Channel  9  r.m.s.  noise  value  of  0. 15- 
0.16  K,  a  value  similar  to  the  0.16K  NEAT  measured  pre¬ 
launch  for  Channel  9  of  the  AMSU-A  instrument  on  Aqua 
(Lambrigtsen,  2003).  Hence  we  will  use  0.16  K  as  our  Chan¬ 
nel  9  noise  floor. 

While  extensive  averaging  of  radiance  variances  can  yield 
reliable  geophysical  variance  estimates  near  or  even  below 
nominal  noise  variances  (Jiang  et  ah,  2004;  Wu,  2004),  here 
we  will  assume  measurement  of  a  single  wave  during  a  sin¬ 
gle  satellite  overpass.  A  conservative  detectability  criterion 
for  such  an  observation  is  a  signal-to-noise  ratio  of  >2,  and 
thus  brightness  temperature  fluctuations  of  ±0.3  K  or  greater. 
While  r.m.s.  temperature  amplitudes  for  gravity  waves  in 


the  extratropical  lower  stratosphere  are  ~1-3K  (e.g.,  Eck¬ 
ermann  et  ah,  1995;  Tsuda  et  ah,  2000),  specific  wave  events 
can  have  peak  amplitudes  as  large  as  10  K  (e.g.,  Hertzog  et 
ah,  2002;  Eckermann  et  ah,  2006a).  Choosing  a  peak  gravity 
wave  temperature  amplitude  7pcak=3  K  at  the  low  end  of  this 
1-10  K  range,  we  require  a  minimum  spectral  visibility  of 
0.1  to  detect  this  wave  at  our  ±0.3  K  lower  limit  for  bright¬ 
ness  temperature.  This  nominal  threshold  for  detectability  is 
highlighted  as  the  bold  contour  in  each  panel  of  Fig.  8. 

Figure  9  plots  the  corresponding  visibility  spectra  for  the 
AMSU-A  on  EOS  Aqua.  The  smaller  footprints  in  Fig.  6b 
yield  Wj(ky,kz)  spectra  that  are  horizontally  elongated 
compared  to  the  corresponding  NOAA  AMSU-A  spectra  in 
Fig.  8,  highlighting  the  Aqua  instmment’s  greater  sensitivity 
to  cross-track  wavelengths  Xy. 

These  2-D  visibility  spectra  at  the  mid-range  scan  angles 
show  small  but  noticeable  asymmetries  about  the  ky—ky—  0 
axes,  which  indicate  via  Eq.  (13)  cross-track  asymmetries  in 
the  response  to  gravity  waves.  These  effects  are  investigated 
further  in  Sect.  4.3. 

4.2  Spectral  response  versus  beam  position  for  a  given 
wave 

Next  we  use  these  Wj(ky,  kz)  spectra  to  investigate  how 
the  AMSU-A  visibility  to  a  2-D  temperature  oscillation  of  a 
given  Xy  and  Xz  varies  at  each  beam  position  j .  Figure  10a 
plots  visibilities  estimated  from  numerical  model  spectra  at 
ky=400  km  and  Xy- 1 2  km.  We  set  ky  positive  and  kz  nega¬ 
tive,  consistent  with  a  gravity  wave  with  an  upward  vertical 
group  velocity  propagating  along  the  positive  Y  axis  direc¬ 
tion  in  Fig.  1. 

We  focus  first  on  visibility  variations  for  the  negative  scan 
angles  f>j  in  Fig.  10a  (j— 1  . . .  15).  For  the  near-zenith  beam 
(7=15),  the  NOAA  AMSU-A  visibility  is  ~  13%,  then  in¬ 
creases  with  increasing  \f}j\  to  a  peak  of  ~13.7%  at  7~ 6- 
8.  From  /'= 6  to  j= 1,  however,  the  response  decreases  with 
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cross-track  wavelength  XY  (km)  cross-track  wavelength  XY  (km) 

-50  km -66.67  km- 100  km  -200  km  200  km  100  km  66.67  km  ±50  km  -66.67  km- 100  km  -200  km  200  km  100  km  66.67  km  50  km 


Fig.  8.  Fourier  Transforms  W;  (ky ,  kz)  computed  from  the  model-generated  2-D  Channel  9  weighting  functions  Wj(Y ,  Z)  for  the  AMSU-A 
on  the  NOAA  15-18  satellites  in  Fig.  4  and  normalized  to  a  peak  response  of  unity,  for  (a)  j=  2  (£2=— 45°),  (b)  j= 6  (/Sg=— 31.67°),  (c)  j=9 
f>S9=— 21.67°9,  and  (d)  j=l5  (/S15=- 1.67°). 


cross-track  wavelength  (km)  cross-track  wavelength  /.v  fkm) 

-50  km -66.67  km- 100  km -200  km  200  km  100  km  66.67  km  ±50  km  -66.67  km- 100  km  -200  km  200  km  100  km  66.67  km  50  km 


Fig.  9.  As  for  Fig.  8,  but  showing  Channel  9  model  2-D  weighting  function  Fourier  Transforms  Wj  ( ky ,  kz)  for  the  AMSU-A  on  EOS  Aqua. 


increasing  \fij\.  This  decreasing  trend  in  visibility  is  due  to 
the  greater  influence  of  cross-track  footprint  widths  (fy)j, 
which  grow  to  diameters  0(100  km)  at  the  outermost  scan 
angles  (see  Fig.  6a)  and  thus  significantly  degrade  detection 
of  this  /,g=400  km  wave.  Note  that  the  corresponding  curve 
for  the  AMSU-A  on  EOS  Aqua  (dotted  curve  in  Fig.  10a) 


shows  larger  visibilities  at  all  j  and  a  peak  visibility  that 
occurs  at  a  larger  \fij\  than  for  the  NOAA  AMSU-A  curve. 
Both  features  are  consistent  with  the  smaller  cross-track  foot¬ 
print  sizes  of  the  Aqua  AMSU-A  at  each  beam  position  (see 
Fig.  6). 
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conjugate  beam  pairs  j,j. 


Fig.  10.  (a)  Spectral  visibilities  W j{ky ,  kg)  as  a  function  of  beam 
position  j  for  ky=400  km  (ky>0)  and  kz=12km  (kz< 0)  for  the 
AMSU-A  on  NOAA  (solid  curve)  and  EOS  Aqua  (dotted  curve), 
(b)  Ratio  RjjJky,  kz )  of  the  negative  scan  angle  response  to  the 
positive  scan  angle  response  in  (a)  for  each  of  the  15  conjugate  beam 
pairs  j  and  /*. 


4.3  Cross-track  asymmetries  in  visibility 

Visibility  curves  at  the  positive  scan  angles  (j  =  16  . . .  30)  in 
Fig.  10a  differ  in  both  magnitude  and  shape  to  those  just  dis¬ 
cussed  at  the  negative  scan  angles  (j— 1  . . .  15).  Figure  10b 
plots  the  ratio  of  these  visibilities  for  each  of  the  15  conju¬ 
gate  beam  pairs  (j,  /'*)•  defined  by  their  identical  off-nadir 
angles  \ftj\=\ftj*\.  From  Eq.  (13),  this  ratio  is  given  by 


Rj,U&Y,kz) 


Wj(ky,  kz) 
Wjjky, kz) 


Wj(ky,kz) 
Wj(-ky, kz) 


(14) 


Asymmetry  is  introduced  here  by  our  wave  structure  and 
our  sign  choices  for  its  wavenumbers  (ky> 0,  kz< 0).  To  first 
order,  a  gravity  wave  is  more  visible  to  a  scanning  instru¬ 
ment  when  its  phase  lines  roughly  coalign  with  the  tilt  angle 
ftj  of  the  line-of-sight  ray  from  the  measurement  point  to 
the  satellite  (e.g.,  Wu  et  ah,  2006).  It  is  straightforward  to 
show  (see  Appendix  A)  that  ft  j=ftj+<f>(ftj):  these  values  are 
plotted  as  the  dotted  curve  in  Fig.  2.  The  phase  lines  of  the 
wave  considered  in  Fig.  10  are  aligned  at  an  off-nadir  angle 
of  arctan(kz/ky)=— 88.28°.  Thus  this  wave  should  be  more 
visible  at  the  negative  scan  angles  than  at  the  positive  scan 
angles,  consistent  with  visibility  ratios  R/jt(ky,  kz)  greater 
than  unity  in  Fig.  10b.  This  also  explains  (to  first  order)  the 
increase  in  these  ratios  with  increasing  |  ftj  |  over  the  y  ~  1 5—5 
range  in  Fig.  10b. 

Interestingly,  however,  this  increase  in  asymmetry  with 
scan  angle  abates  at  y  ~ 5  and  the  cross-track  asymmetry  then 
decreases  with  increasing  \ftj  at  the  outermost  scan  angles 
( j—  1  . .  .5).  This  feature  is  somewhat  surprising  given  that 
the  outermost  beam  position  pairs  (  j,  y*)  have  the  great¬ 
est  differences  in  their  respective  scan  angles  ftj  and  ftjx , 
and  thus  might  be  expected  to  exhibit  the  largest  cross-track 
asymmetries  in  their  gravity  wave  visibilities.  That  they  do 
not  merits  analysis. 


In  studying  large  cross-track  asymmetries  in  the  response 
of  MLS  saturated  limb-track  radiances  to  gravity  waves, 
McLandress  et  al.  (2000)  and  Jiang  et  al.  (2004)  introduced 
the  concept  of  an  effective  line-of-sight  angle,  ft.  The  com¬ 
bined  effects  of  vertical  variations  in  absorption  and  antenna 
spreading  across  the  MLS  viewing  direction  combined  to 
yield  a  final  2-D  MLS  weighting  function  that  was  not  sym¬ 
metric  about  the  line  of  sight  direction  ft,  but  was  instead 
approximately  symmetric  about  an  effective  line  of  sight  di¬ 
rection  ft  that  was  more  horizontally  aligned  than  ft  (Wu  and 
Waters,  1996a).  Gaussian  analytical  approximations  to  the 
MLS  vertical  weighting  functions  and  antenna  pattern  en¬ 
abled  McLandress  et  al.  (2000)  to  evaluate  ft  analytically. 
Rotating  the  (Y,Z)  axes  to  new  tilted  axes  (Y'.Z'),  such  that 
Y'  was  aligned  along  this  MLS  effective  line-of-sight  direc¬ 
tion  ft,  led  to  a  separable  2-D  MLS  weighting  function  with 
elliptical  contours  whose  long  axis  was  aligned  along  Y'  and 
whose  short  axis  was  aligned  along  Z'.  Peak  responses  oc¬ 
curred  for  gravity  waves  whose  phase  lines  were  aligned  par¬ 
allel  to  Y\  since  their  oscillations  (orthogonal  to  phase  lines 
along  their  wavenumber  vector)  occur  along  the  shortest  Z' 
axis  of  the  weighting  function  and  thus  suffer  the  least  smear¬ 
ing.  Hence  the  term  “effective  line  of  sight”  for  this  MLS  Y' 
axis. 

We  can  derive  similar  angles  ftj  for  the  AMSU-A  weight¬ 
ing  functions,  though  the  nomenclature  “effective  line  of 
sight”  is  not  apt  here,  since,  unlike  limb  viewers  like  MLS, 
ft  j  and  ft  j  are  very  different  for  sub  limb  viewers  like  AMSU- 
A.  Hence  we  refer  to  ftj  as  the  cross-track  tilt  angle  of  the 
AMSU-A  weighting  function,  and  Y’-  as  its  cross-track  tilt 
axis.  Also,  unlike  the  MLS  case,  here  we  cannot  evaluate 
ftj  analytically  since  our  AMSU-A  weighting  functions  are 
numerical,  nonGaussian  in  the  vertical  and  somewhat  asym¬ 
metric  cross-track  due  to  spherical  effects.  Nonetheless,  we 
can  estimate  ftj  numerically  as  follows. 

First,  we  note  from  Jiang  et  al.  (2004)  that  their  MLS  ef¬ 
fective  line-of-sight  axis  Y'  can  be  defined  approximately  as 
the  line  interconnecting  local  maxima  in  vertical  profiles  of 
the  2-D  weighting  functions  at  different  cross  track  positions 
Y :  see,  e.g.,  their  Fig.  8c.  Thus,  we  estimate  cross-track  tilt 
angles  ftj  by  evaluating  altitudes  Z  of  the  local  maxima  in 
vertical  profiles  of  W  j  (Y .  Z)  at  a  series  of  different  Y  values. 
The  resulting  locus  of  (T,  Z)  coordinates  specifies  the  cross¬ 
track  tilt  axis  Yj. :  the  result  of  such  a  calculation  for  the  ;=7 
Aqua  AMSU-A  2-D  weighting  function  is  plotted  as  the  gray 
thick  solid  line  in  Fig.  11a.  The  off-nadir  tilt  angle  of  this  line 
is  ftj.  One  can  see  from  Fig.  11a  that  the  width  of  the  2-D 
weighting  functions  along  a  Z'-  axis  that  lies  orthogonal  to 
this  solid  gray  Y'.  axis  is  narrower  than  the  width  along  the 
dashed  line-of-sight  axis,  for  example.  Thus,  waves  whose 
phase  lines  align  parallel  to  Y'.  should  be  more  visible  than 
those  whose  phases  are  aligned  differently. 

Figure  lib  plots  our  numerical  estimates  of  tilt  angles  ftj 
as  a  function  of  beam  position  j  for  both  the  NOAA  and 
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(a)  EOS  Aqua  j=  7 


cross-track  distance  Y  (km) 


3029282726252423222120191817161514131211109  87654321 

beam  number  j 


Fig.  11.  (a)  Wj  (T,  Z)  at  j=l  for  AMSU-A  on  EOS  Aqua.  Contours 
are  values  relative  to  the  peak  response,  expressed  as  a  percentage. 
The  50%  contour  is  shown  in  bold.  Dotted  gray  line  shows  the  line- 
of-sight  ray  to  the  satellite.  Solid  gray  curve  connects  local  maxima 
in  vertical  cross  sections  through  the  weighting  function  at  different 
Y,  and  specifies  the  cross-track  tilt  axis  Y'-.  The  off-nadir  tilt  angle 

of  this  gray  line  is  fi : .  (b)  Pj  angles  computed  as  in  (a)  for  all  the 
weighting  functions  Wj  ( Y ,  Z )  for  both  the  NOAA  (solid  curve)  and 
EOS  Aqua  (dotted  curve)  AMSU-A. 


Aqua  AMSU-A  weighting  functions.  These  angles  are  all 
close  to  —90°  (i.e.  very  little  axis  tilt).  On  progressing  out¬ 
wards  from  nadir  in  Fig.  lib,  these  angles  become  progres¬ 
sively  more  tilted  out  to  the  mid-range  scan  angles,  then  level 
off  and  start  to  return  towards  more  horizontal  Y' ■  axes  at  the 
outermost  scan  angles.  This  return  to  symmetry  at  the  out¬ 
ermost  scan  angles  is  due  mostly  to  much  wider  horizontal 
footprint  diameters  (fy)  j  from  antenna  spreading.  For  ex¬ 
ample,  the  analytical  expressions  of  McFandress  et  al.  (2000) 
and  Jiang  et  al.  (2004)  show  that  as  the  ratio  of  the  widths 
of  the  antenna  spreading  to  vertical  absorption  (their  y~1^2 
parameter)  increases  at  the  measurement  point,  the  MFS  ef¬ 
fective  line-of-sight  axis  Y'  gets  pushed  further  towards  the 


40  20  ^ 

grees) 

J  -20  -40 

(a)  :  >> 

f— ^  \  '  : 

A,y=200  km 

km 

\\  - 

\\ 

j 1 

\ 

\ 

V 

3028262422201816141210  8  6  4  2 
beam  number  j 


0  10  40 


conjugate  beam  pairs  j ,  j. 


Fig.  12.  As  for  Fig.  10,  but  for  Xy=200km  and  kz=25  km. 

horizontal.  This  width  ratio  becomes  large  at  the  outermost 
AMSU-A  scan  angles  since  antenna  spreading  widths  in¬ 
crease  significantly  here  due  to  longer  ray  path  lengths  from 
the  satellite  to  the  measurement  point.  Thus  the  trend  back  to 
symmetric  responses  at  far  off-nadir  scan  angles  in  Fig.  10b 
is  consistent  with  the  transition  to  less-tilted  Y ' ■  axes  for  the 
weighting  functions  in  Fig.  1  lb. 

Figure  12  shows  visibility  results  for  a  different  wave  of 
/.)  =200  km  and  a 7 =25  km,  with  ky> 0  and  kz< 0  as  before. 
Peak  visibilities  in  Fig.  12a  are  larger  (~40-45%)  than  those 
in  Fig.  10a  due  to  the  longer  vertical  wavelength  which  is 
less  smeared  by  the  vertical  weighting  functions.  Since  the 
shorter  cross-track  wavelength  is  more  prone  to  degraded 
detection  by  the  broader  off-nadir  footprints,  however,  vis¬ 
ibilities  start  decreasing  with  increasing  fi  j  j  earlier  than  in 
Fig.  10a.  As  in  Fig.  10b,  the  cross-track  visibility  ratios 
RjJt  (kv ,  kz)  in  Fig.  12b  increase  with  increasing  be¬ 
fore  rolling  off  towards  more  symmetric  ratios  at  the  outer¬ 
most  scan  angles. 


5  Three-dimensional  forward  model  simulations 

Here  we  simulate  gravity  wave  detection  more  directly  and 
in  three  dimensions  by  numerically  evaluating  Eq.  (11)  in 
the  presence  of  a  3-D  gravity  wave  temperature  oscillation 
using  the  3-D  weighting  functions  from  our  AMSU-A  radi¬ 
ance  acquisition  model.  Our  goal  is  to  assess  whether  the 
two-dimensional  horizontal  coverage  of  the  AMSU-A  mea¬ 
surements  in  Fig.  6  can  “image”  the  resolved  gravity  wave 
radiance  perturbations  in  the  horizontal. 

We  begin  by  specifying  an  idealized  infinite  three- 
dimensional  train  of  monochromatic  gravity  wave  tempera¬ 
ture  oscillations  of  the  form 

T\X,  T,  Z)  =  Tpeak (Z)  cos  (kxX  +  kyY  +  kzZ) ,  (15) 

where  7pCak  is  the  peak  wave  temperature  amplitude  (po¬ 
tentially  height  dependent).  We  specify  the  wave’s  vertical 
wavelength  Xz—2it /\kz\.  horizontal  wavelength  k/,=27r/k/i, 
and  horizontal  propagation  azimuth  (p  with  respect  to  the 
AMSU-A  viewing  geometry  axes  (X.Y.Z)  in  Fig.  1,  such 
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that  kx—^2jt/Xx—kh  cos  (p  and  ky=±2n  /Xy=kh  sin<£>.  We 
then  sample  this  temperature  structure  from  orbit  using  the 
cyclical  AMSU-A  scanning  pattern,  acquiring  brightness 
temperature  perturbations  T'B (Xj ,  Yj )  by  evaluating  Eq.  (11) 
using  Eq.  (15)  as  the  3-D  temperature  field.  This  calculation 
assumes  that  the  background  brightness  temperature  upon 
which  these  perturbations  are  superimposed  (e.g..  Fig.  7)  can 
be  perfectly  extracted  from  the  measurements:  see  Ecker¬ 
mann  et  al.  (2006b)  for  examples  of  how  this  is  done  (imper¬ 
fectly)  in  practice.  Equation  (1 1)  is  evaluated  by  rectangular 
integration  after  interpolating  the  model-generated  weight¬ 
ing  functions  Wj(X,  Y,  Z)  onto  a  regular  ( X ,  Y,  Z )  grid  with 
10  km  horizontal  and  0.5  km  vertical  grid  spacings. 

Figure  13  shows  NOAA  AMSU-A  results  for  a  wave  of 
A/, =400  km  and  ),y=  1 2  km,  aligned  at  three  different  hori¬ 
zontal  propagation  azimuths  q>.  A  constant  peak  amplitude 
7peak=5  K  was  used.  Since  there  is  no  explicit  dependence 
of  these  results  on  wave  amplitude,  instead  of  plotting  ab¬ 
solute  brightness  temperatures  T'B(Xj,  Y /)  we  instead  plot 
normalized  visibility  perturbations  T'b(Xj ,  T/)/7peak,  to  fa¬ 
cilitate  direct  comparisons  with  the  2-D  spectral  visibilities 
in  Sect.  4. 

The  2-D  weighting  function  spectra  for  this  wave  at 
i p—90°  were  profiled  in  Fig.  10a.  Figure  13a  shows  re¬ 
sults  from  our  3-D  sampling  of  this  wave  for  <£>=80°.  As  in 
Fig.  10b,  there  is  weak  but  distinct  cross-track  asymmetry  in 
the  visibilities  about  the  subsatellite  point,  with  radiance  vis¬ 
ibility  perturbation  amplitudes  of  ~  1 3%,  a  peak  response  at 
T~— 500  km  0'~6),  then  a  falloff  to  peak  visibilities  nearer 
10%  at  the  outermost  scan  angles.  For  our  5  K  peak  wave 
amplitude,  visibility  oscillations  of  ±13%  imply  brightness 
temperature  perturbations  T'B  of  ±0.65  K,  well  above  our 
conservative  nominal  detectability  threshold  of  ±0.3  K  dis¬ 
cussed  in  Sect.  4.1.  Similar  amplitudes  occur  for  the  other 
two  propagation  directions  q>  in  Figs.  10b  and  c.  The  model 
results  clearly  indicate  that  NOAA  AMSU-A  Channel  9  ra¬ 
diances  can  image  this  particular  gravity  wave,  regardless  of 
its  propagation  angle  with  respect  to  the  scanning  geome¬ 
try  during  the  satellite  overpass.  The  same  holds  for  EOS 
Aqua,  since  its  visibilities  are  greater  at  all  beam  positions 
than  those  of  NOAA  AMSU-A  (see  Fig.  10a). 

However,  the  imaged  wave  structure  in  Fig.  13  changes 
noticeably  with  varying  propagation  direction  <£>.  When 
the  gravity  wave  horizontal  wavenumber  Kh—(kx,ky)  is 
aligned  roughly  parallel  to  the  cross-track  Y  axis  (<£>=80°: 
Fig.  13a),  horizontal  wave  phase  lines  (orthogonal  to  Kh)  are 
well  imaged  by  the  near-nadir  scans,  but  become  poorly  im¬ 
aged  at  the  outer  scan  angles  due  to  the  larger  cross-track 
footprint  diameters  (fy)j-  This  is  entirely  consistent  with 
the  falloff  in  2-D  spectral  visibility  for  this  wave  (^=80°) 
at  large  \fij  j  noted  in  Fig.  10a.  For  ^>=45°  (Fig.  13b)  and 
<p=350°  (Fig.  13c),  the  horizontal  wave  structure  is  imaged 
better  at  the  outermost  scan  angles  than  for  ^=80°,  since 
these  wave  orientations  project  longer  wavelength  compo¬ 
nents  cross-track  which  are  more  easily  resolved  by  the  rel¬ 


atively  wider  cross-track  footprint  diameters  (fy)j-  From 
Fig.  5,  the  along-track  footprints  (fx)j  at  the  far  off-nadir 
beam  positions  are  ~56%  the  size  of  the  cross-track  foot¬ 
prints  ( fy)j ,  and  are  oversampled  in  the  along-track  ( X)  di¬ 
rection  (see  Fig.  6a).  Both  properties  allow  shorter  horizon¬ 
tal  wavelength  structure  to  be  imaged  along-track  compared 
to  cross-track  at  the  outermost  scan  angles,  leading  to  better 
imaging  of  these  latter  wave  orientations  at  the  far  off-nadir 
locations  in  Figs.  13b  and  c  than  in  Fig.  13a. 

The  imaged  wave  fields  yield  phase  and  wavelength  infor¬ 
mation.  The  imaged  gravity  wave  phase  lines  in  Figs.  13b 
and  c  show  a  noticeable  curving  or  bowing,  away  from  their 
intrinsic  linear  forms.  This  is  produced  by  the  upward  dis¬ 
placement  in  the  peaks  of  the  vertical  weighting  functions 
with  increasing  scan  angle  in  Figs.  3  and  4,  due  to  the  limb 
effect.  This  means  that  wave  phase  sampled  at  different 
beam  positions  j  varies  not  just  due  to  the  change  in  hor¬ 
izontal  viewing  location  (Xj,  Yj),  but  also  with  the  verti¬ 
cal  change  in  the  altitude  of  the  weighting  function  peaks. 
Since  the  former  effect  is  dominant,  the  radiance  perturba¬ 
tions  do  reproduce  an  approximate  horizontal  cross  section 
of  the  wave  oscillation  pattern,  such  that  the  wavelengths  in 
Fig.  13  reproduce  the  actual  400km  horizontal  wavelength  of 
this  wave  quite  well.  The  curvature  added  to  these  phase  lines 
is  the  smaller  secondary  influence  of  vertical  phase  changes 
induced  by  slightly  different  peak  heights  for  the  weighting 
functions  W/(X,  Y,  Z)  at  each  beam  position  j.  Thus  this 
curving  of  phase  lines  is  the  corresponding  “limb  effect”  dis¬ 
tortion  for  AMSU-A  gravity  wave  measurements  that  “limb 
darkening”  and  “limb  brightening”  (Kidder  and  Vonder  Haar, 
1995)  are  for  AMSU-A  background  temperature  measure¬ 
ments  (see  Figs.  7b  and  c). 

Figure  14  shows  results  for  the  same  three  propagation  az¬ 
imuths  as  in  Fig.  13,  but  for  a  gravity  wave  with  a  shorter  hor¬ 
izontal  wavelength  (A/,  =200  km)  and  longer  vertical  wave¬ 
length  (Az=25km).  2-D  NOAA  AMSU-A  visibilities  for 
this  wave  (^>=90°)  in  Fig.  12a  showed  peak  visibility  am¬ 
plitudes  of  ~43%,  decreasing  rapidly  with  increasing  \fi  ;  |  at 
the  outermost  scan  angles  to  values  of  ~10-20%.  The  <£=80° 
forward-modeled  swath  imagery  for  this  wave  in  Fig.  14a  re¬ 
produces  these  same  trends  in  visibility. 

Like  the  A/, =400  km  wave  in  Fig.  13,  the  horizontal  imag¬ 
ing  of  this  wave’s  oscillations  degrades  with  increasing  |/!y| 
more  seriously  for  <£=80°  (Fig.  14a)  than  for  the  other  prop¬ 
agation  azimuths  in  Figs.  14b  and  c.  When  this  A/,=200km 
wave  propagates  cross-track  (Fig.  14a),  there  is  a  highly 
nonuniform  amplitude  response  cross-track  in  the  imaged 
wave  structure  since  only  the  narrow  near-nadir  AMSU-A 
footprints  resolve  it  well.  Conversely,  when  this  wave  prop¬ 
agates  more  nearly  along  track  (Fig.  14c),  the  narrower  foot¬ 
prints  and  oversampling  along-track  at  the  outermost  scan 
angles  yield  a  more  uniform  imaged  response  to  this  wave  in 
the  total  swath  coverage.  As  in  Fig.  13,  we  also  see  curving 
of  imaged  wave  phase  lines  for  this  wave,  most  noticeably  at 
^=45°  in  Fig.  14b. 
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(a)  A,„=400km,  3.z=12km,  (p=80°  (b)  3.h=400km,  X2=12km,  (p=45°  (c)  ?ih=400km,  3.z=12km,  cp=350°  Tg/T,,^ 
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Fig.  13.  (a)  Relative  brightness  temperature  perturbations  (visibilities)  T^(Xj,  Y j)/T^ea ^  resulting  from  forward-modeled  NOAA  AMSU-A 
Channel  9  sampling  of  a  gravity  wave  in  Eq.  (15)  with  Tpeak=5  K,  k/,=400km,  >-2=12  km,  and  tp  values  of  (a)  80°,  (b)  45°,  and  (c)  350°. 
The  white  vector  at  the  center  of  each  plot  shows  this  direction  of  horizontal  wave  propagation  <p.  The  color  scale  is  the  visibility  amplitude 
expressed  as  a  percentage,  indicating  a  peak  sensitivity  of  ~13%,  or  a  maximum  brightness  temperature  perturbation  of  ~0.65  K. 


(a)  /.ll=200krn,  /,z=25km,  (p=80°  (b)  /.h=200km.  /v=25km.  (p=45°  (c)  /-=200km.  /v=25km,  cp=350° 
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Fig.  14.  Same  presentation  as  in  Fig.  13,  but  for  a  gravity  wave  with  k/,  =200 km,  kz=25  km. 


6  Discussion 

To  investigate  the  sensitivity  of  AMSU-A  to  atmospheric 
gravity  waves,  we  have  developed  a  simple  numerical  model 
of  its  in-orbit  radiance  acquisition.  Of  the  six  available 
stratospheric  temperature  channels,  we  focused  on  Chan¬ 
nel  9.  Since  Channel  9  views  a  single  O2  line  and  its  1- 
D  vertical  weighting  functions  peak  at  ~60-90hPa,  its  ra¬ 
diative  transfer  should  more  closely  approximate  the  purely 
pressure-broadened  single  Lorentz-line  absorption  model  we 
developed  in  Sect.  3.  Furthermore,  Channel  9  has  slightly 
narrower  vertical  weighting  functions  than  the  other  strato¬ 
spheric  channels  (Goldberg  et  ah,  2001),  making  it  (theo¬ 
retically)  more  sensitive  to  gravity  waves.  Its  lower  altitude 
coverage  also  enables  us  to  model  any  waves  seen  in  these 
channel  radiances  with  global  and  regional  numerical  pre¬ 
diction  models,  which  can  simulate  both  wave  generation  in 
the  troposphere  and  the  propagation  of  these  waves  into  the 
lower  stratosphere  (e.g.,  Wu  and  Zhang,  2004;  Eckermann  et 
ah,  2006b). 

Using  both  spectral  analysis  of  our  model-derived  2-D 
weighting  functions  and  numerical  forward  modeling  using 


our  3-D  weighting  functions  applied  to  3-D  gravity  wave 
fields,  we  have  shown  that  certain  types  of  gravity  wave  os¬ 
cillations  appear  in  Channel  9  radiances  as  brightness  tem¬ 
perature  fluctuations  with  amplitudes  well  in  excess  of  antic¬ 
ipated  noise  floors.  While  our  results  here  are  specific  to  this 
single  channel,  they  argue  for  gravity  wave  detection  in  radi¬ 
ances  from  the  other  higher-altitude  AMSU-A  Channels  10- 
14  as  well,  given  the  identical  scanning  pattern  and  similar 
footprint  diameters  and  vertical  weighting  function  widths 
(Goldberg  et  ah,  2001).  While  noise  floors  (NEAT)  increase 
somewhat  for  the  higher-altitude  channels  (e.g.,  Lambrigt- 
sen,  2003),  so  too  do  typical  gravity  wave  temperature  am¬ 
plitudes  and  thus  anticipated  brightness  temperature  pertur¬ 
bations. 

These  findings  support  the  observational  results  of  Wu 
(2004),  who  analyzed  radiance  perturbations  in  AMSU-A 
Channel  13  radiances  (which  peak  at  ~5  hPa)  and  found  en¬ 
hancements  in  radiance  variances  that  had  very  similar  geo¬ 
graphical  distributions  to  gravity  wave-related  enhancements 
in  MLS  radiances  for  the  same  observation  period.  In  aver¬ 
aging  AMSU-A  radiances,  Wu  (2004)  binned  them  within 
6  different  groups  classified  according  to  beam  position. 
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cross-track  wavenumber  IAryl  (eye  km'1) 


Fig.  15.  Black  contours  show  the  spectral  variance  visibility  ratio 

R 7  •  iky  >  kz)  as  a  function  of  cross-track  and  vertical  wavenumber 
J’J* 

for  the  NOAA  AMSU-A  j=  4,  j*=21  conjugate  beam  pairing.  Gray 
lines  overlay  absolute  visibility  contours  of  Wj  (ky ,  ky)  for  the  j=  4 
beam. 


Groups  1  and  6  used  data  from  the  10  farthest  off-nadir  beam 
positions  (/=  1-5  and  j  =26-30),  with  group  1  incorporating 
data  in  which  the  beam  pointed  east  (Egl),  and  group  6  using 
the  remaining  data  in  which  the  beam  pointed  west  (Wg6). 
We  denote  these  brightness  temperature  variances  <7^gl  and 

°Wg6’  respectively.  On  mapping  these  variances  in  the  winter 
extratropical  Southern  Hemisphere,  Wu  (2004)  found  signif¬ 
icant  differences  over  the  Andes  and  New  Zealand,  with  crgal 
typically  a  factor  of  2  larger  than  er^,g6.  This  suggests  cross¬ 
track  asymmetries  in  the  response  to  the  stratospheric  gravity 
waves  over  these  land  masses. 

In  our  analysis  of  modeled  NOAA  AMSU-A  Channel  9 
weighting  functions,  cross-track  asymmetries  in  wave  am¬ 
plitude  on  the  order  of  10-15%  (20-30%  in  variance)  were 
noted  for  two  particular  wave  case  studies.  To  study  asym¬ 
metries  over  a  wider  range  of  wavelength  pairings.  Fig.  15 
plots  cross-track  variance  ratios  R 7  .  (ky,  ky)  for  the  NOAA 
AMSU-A  y=4,  y*=27  conjugate  beam  pairing.  Within  re¬ 
gions  of  >10%  absolute  visibility  for  7  =4  (the  thick  solid 
gray  contour  in  Fig.  15),  we  see  cross-track  variance  ratios 
of  up  to  1.5,  below  the  factor  of  2  reported  by  Wu  (2004) 
for  Channel  13.  Factors  of  2  do  occur  in  Fig.  15,  but  at 
larger  wavenumbers  where  the  visibilities  are  so  small  that 
we  would  anticipate  little  or  no  radiance  signal  above  the 
noise  floor  unless  those  waves  had  very  large  amplitudes 
7 peak  •  The  larger  anisotropies  reported  by  Wu  (2004)  may  re¬ 
sult  from  his  cross-track  fitting  method  used  to  isolated  wave 
fluctuations,  which  tends  to  isolate  waves  propagating  cross¬ 
track  better  than  those  propagating  along  track,  since  the  lat¬ 
ter  project  longer  wavelengths  cross-track  (see  Fig.  13c)  that 
are  harder  to  separate  from  mean  cross-track  limb  effects  (see 
Fig.  7c).  In  addition,  differences  in  the  weighting  functions 
between  Channels  13  and  9,  or  shortcomings  in  our  simpli¬ 


fied  treatment  of  radiative  transfer  and  antenna  specifics,  may 
also  yield  underestimates  in  cross-track  anisotropies. 

Our  3-D  forward  model  simulations  in  Sect.  5  also  qualita¬ 
tively  confirm  preliminary  findings  of  Wu  and  Zhang  (2004) 
that  the  two-dimensional  horizontal  coverage  provided  by 
AMSU-A’s  cross-track  scanning  pattern  yields  horizontal 
images  of  resolved  gravity  wave  oscillations.  A  mesoscale 
model  simulation  by  Wu  and  Zhang  (2004)  on  the  days  and 
at  the  locations  of  their  wave  measurements  produced  gravity 
waves  in  the  lower  stratosphere  with  horizontal  wavelengths 
of  ~300-500  km  and  vertical  wavelengths  of  ~7 — 15  km  (see 
their  Fig.  13).  The  approximate  midpoint  of  this  range  corre¬ 
sponds  to  the  /,/,=400  km  and  /.y= 12  km  wave  we  considered 
in  both  2-D  and  3-D  model  simulations  (Figs.  10  and  13,  re¬ 
spectively).  Unfortunately,  Wu  and  Zhang  (2004)  only  plot¬ 
ted  wave-induced  divergence  perturbations  from  their  model 
runs,  but  they  quote  the  largest  temperature  amplitudes  in 
their  model  to  be  ~5  K.  Our  earlier  modeling  predicted  peak 
visibilities  of  ~  1 3%  for  this  wave  and  a  maximum  Channel  9 
brightness  temperature  amplitude  of  ~0.65  K  for  a  5  K  wave 
temperature  amplitude.  This  0.65  K  value  is  somewhat  larger 
than  the  ±0. 3-0.5  K  Channel  9  brightness  temperature  oscil¬ 
lations  reported  by  Wu  and  Zhang  (2004)  (see  their  Figs.  5, 
7  and  10).  However,  given  the  large  uncertainties  in  the 
actual  temperature  amplitudes  and  vertical  wavelengths  for 
these  wave  oscillations  and  the  strong  senstivities  of  wave- 
induced  Channel  9  radiance  oscillations  to  uncertainties  in 
both  of  these  parameters  (e.g..  Fig.  8),  these  differences  can¬ 
not  be  considered  definitive.  Definitive  assessments  require 
direct  forward  modeling  of  the  full  3-D  wave  temperature 
fields  into  corresponding  model-specified  Channel  9  bright¬ 
ness  temperatures,  whose  perturbation  structures  can  then 
be  compared  directly  to  those  observed.  Such  comparisons 
are  the  focus  of  the  companion  paper  of  Eckermann  et  al. 
(2006b).  For  a  large-amplitude  stratospheric  gravity  wave 
over  southern  Scandinavia,  they  report  close  agreement  be¬ 
tween  observed  Channel  9  radiance  perturbations  and  those 
forward  modeled  using  our  3-D  weighting  functions  and  3-D 
temperature  oscillations  from  numerical  weather  prediction 
models. 

Our  model  results  for  AMSU-A  also  provide  useful  guid¬ 
ance  for  assessing  the  ability  of  future  scanning  microwave 
sensors  to  detect  gravity  waves.  For  example,  the  Ad¬ 
vanced  Technology  Microwave  Sounder  (ATMS),  slated  to 
fly  first  on  the  National  Polar-orbiting  Operational  Envi¬ 
ronmental  Satellite  System  (NPOESS)  Preparatory  Project 
(NPP)  and  then  on  operational  NPOESS  platforms,  will  have 
similar  scanning  and  temperature  channel  sensitivities  to 
AMSU.  Its  smaller  nominal  nadir  surface  footprint  diame¬ 
ters  of  ~33  km  should  make  ATMS  more  sensitive  to  shorter 
horizontal  wavelengths.  The  Special  Sensor  Microwave  Im¬ 
ager/Sounder  (SSMIS)  on  the  Defense  Meteorological  Satel¬ 
lite  Program  (DMSP)  satellite  launched  in  October  2003, 
which  uses  a  conical  scanning  strategy  (see  Rosenkranz 
et  al.,  1997),  yields  radiances  in  its  lower  stratospheric 
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temperature  Channels  6  and  7  with  footprint  diameters  of 
-37.5  km  and  nominal  NEAT-0.4-0.5  K.  Our  AMSU-A 
modeling  indicates  that  the  better  horizontal  resolution  of  SS- 
MIS  coupled  with  its  narrower  vertical  weighting  functions 
due  to  the  fixed  45°  off-nadir  beam  position  should  allow  this 
instrument  to  detect  gravity  wave  fluctuations  in  its  Channel 
6-7  radiances,  despite  somewhat  higher  nominal  noise  floors 
compared  to  AMSU-A.  Its  fixed  scan  angle  might  also  re¬ 
duce  the  limb-effect  distortion  of  phase  lines  noted  here  for 
AMSU-A. 


7  Summary  and  conclusions 

The  modeling  work  outlined  here  has  shown  that  gravity 
waves  with  long  vertical  wavelengths  (kz>10km),  long  hor¬ 
izontal  wavelengths  (k/?>  150-200  km),  and  typical  lower 
stratospheric  temperature  amplitudes  (Tpeak^l-3 K)  are  re¬ 
solved  as  perturbations  to  the  lower  stratospheric  microwave 
radiances  acquired  by  AMSU-A  Channel  9  on  both  the 
NOAA  meteorological  satellites  and  NASA’s  EOS  Aqua 
satellite.  One  entire  AMSU-A  scan  cycle  yields  radiance 
measurements  at  30  cross-track  measurement  locations  span¬ 
ning  ±800-1 100  km  either  side  of  the  satellite  ground  track. 
As  the  scan  cycle  repeats  and  the  satellite  motion  sweeps  out 
long  two-dimensional  “pushbroom”  radiance  images  along 
track,  the  horizontal  structure  of  these  resolved  waves  is  im¬ 
aged  in  these  radiance  maps. 

Our  modeling  has  revealed  some  interesting  instrumen¬ 
tal  effects  that  should  aid  interpretation  of  the  gravity  wave 
signals  in  these  data.  First,  the  amplitude  of  the  gravity 
wave’s  radiance  response  varies  systematically  with  cross¬ 
track  beam  position  j .  As  the  cross-track  viewing  angle  Pj 
changes  from  near-nadir  to  further  off-nadir,  the  AMSU-A 
weighting  functions  become  broader  cross-track  and  slightly 
narrower  vertically  due  to  the  combined  effects  of  antenna 
spreading,  the  limb  effect  and  Earth  curvature.  These  two 
width  changes  combine  to  yield  2-D  weighting  functions  that 
are  tilted  farthest  off-horizontal  at  the  mid-range  scan  an¬ 
gles.  These  tilts  lead  to  asymmetric  responses  to  gravity 
waves  when  viewed  at  conjugate  beam  positions  j  and  /* 
(Pjt——Pj),  with  amplitude  differences  typically  on  the  or¬ 
der  of  10-15%.  The  narrower  footprint  diameters  and  greater 
sampling  rates  along  track  compared  to  cross  track  at  the  out¬ 
ermost  scan  angles  allow  gravity  waves  to  be  imaged  more 
effectively  when  they  propagate  along  track.  The  increase  in 
altitude  of  the  weighting  function  peak  with  increasing  scan 
angle  |/3y  |  introduces  a  distorting  curvature  to  the  gravity 
wave  phase  lines  imaged  in  these  radiances,  which  can  affect 
observational  estimates  of  wave  propagation  directions  or  3- 
D  wave  structures.  Nonetheless,  for  the  waves  we  modeled, 
the  radiance  maps  provide  a  fairly  accurate  (and  valuable) 
direct  measurement  of  the  wave’s  horizontal  wavelength. 

In  a  companion  paper,  Eckermann  et  al.  (2006b)  apply  the 
3-D  weighting  function  results  and  forward  model  methodol¬ 


Fig.  16.  AMSU-A  observational  geometry.  The  instrument  views 
the  atmosphere  at  a  channel  height  Zq  at  a  scan  angle  Pj .  The  line- 
of-sight  distance  from  this  measurement  point  to  the  satellite  is  s 
and  from  the  surface  to  this  measurement  point  is  s.  Y j  is  the  cross¬ 
track  swath  distance  from  the  subsatellite  point  to  the  measurement 
point.  Other  symbols  are  defined  elsewhere  in  the  text. 


ogy  we  have  developed  here  to  an  observational  case  study  of 
a  gravity  wave  over  southern  Scandinavia.  This  study  com¬ 
bines  AMSU-A  Channel  9  radiance  imagery  and  3-D  numer¬ 
ical  model  simulations  of  the  wave  field,  in  an  attempt  to 
fully  characterize  the  wave  and,  via  forward  modeling,  pro¬ 
vides  an  observational  validation  of  this  paper’s  theoretical 
predictions  of  anticipated  AMSU-A  radiance  responses  to 
specific  gravity  wave-induced  3-D  temperature  structure. 


Appendix  A 


AMSU-A  observational  geometry  and  derivation  of 
footprint  diameter  ratios 


The  geometry  is  depicted  in  Fig.  16.  From  the  law  of  sines. 


P  j  —  jt  —  i/r  =  arcsin 


(Re  +  Zsat)  sin  Pj 
Re  +  Zc 


(Al) 


This  is  the  off-nadir  cross-track  angle  that  the  beam  makes  at 
a  measurement  point  at  altitude  Zc,  differing  from  Pj  due 
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to  Earth  curvature.  Since  <j>+Pj+\jr—7T,  then  substituting 
Eq.  (Al)  yields  Eq.  (10),  which  in  turn  implies  that 

Pj  =  Pj  +  4>(Pj)-  (A2) 


The  cross-ttack  distance  Y,  from  the  subsatellite  point 
to  the  measurement  point  at  the  altitude  Zc  is  then  just 
( Re+Zc)4>(Pj ).  Then  we  can  estimate  the  cross-track  half¬ 
power  footprint  diameters  to  be 


(fy)j  —  [Re  +  Zc] 


Phpbw  \ 
2  ) 


/^HPBW  \ 
2  ) 


(A3) 


Along  track,  we  estimate  the  footprint  diameter  using  the 
line-of-sight  ray  distance  s  between  the  satellite  and  mea¬ 
surement  point  and  the  width  of  the  antenna  polar  diagram. 
The  law  of  sines  yields 


(Re  +  Zc)  sin  d>(/L) 

s  =  — - — - 2-^4,  (A4) 

sin  j 

whereupon,  using  a  flat  Earth  approximation  (adequate  here 
given  the  small  antenna  width  /Ihpbw  and  no  along-track  tilt¬ 
ing  of  the  beam). 


(fx)  j  =  2?  tan 


^HPBW 

2 


(A5) 


Equations  (A3)-(A5)  yield  Eq.  (9),  the  along-track  to 
cross-track  footprint  width  ratio. 
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